f = recipient ~ recipient_lag + p_polity2 + gdppc + pop + factor(year)*factor(agency)
f1 = update(f, '.~.+alliance*independence')
f2 = update(f, '.~.+trade*independence')
f3 = update(f, '.~.+infant_mortality*independence')
f = oda ~ oda_lag + p_polity2 + gdppc + pop + factor(year)*factor(agency)
f4 = update(f, '.~.+alliance*independence')
f5 = update(f, '.~.+trade*independence')
f6 = update(f, '.~.+infant_mortality*independence')

mod1 = glm(f1, data=dat, family=binomial())
mod2 = glm(f2, data=dat, family=binomial())
mod3 = glm(f3, data=dat, family=binomial())
mod4 = lm(f4, dat[dat$oda>0,])
mod5 = lm(f5, dat[dat$oda>0,])
mod6 = lm(f6, dat[dat$oda>0,])
selection = list(mod1, mod2, mod3, mod4, mod5, mod6)

print(screenreg(selection, omit.coef='factor', digits=3))
